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We expose at a previously unknown level of detail the features of the dynamics of trajectories that 
either evolve towards the Feigenbaum attractor or are captured by its matching repellor. Amongst 
these features are the following: i) The set of preimages of the attractor and of the repellor are 
embedded (dense) into each other, ii) The preimage layout is obtained as the limiting form of 
the rank structure of the fractal boundaries between attractor and repellor positions for the family 
of supercycle attractors. iii) The joint set of preimages for each case form an infinite number of 
families of well-defined phase-space gaps in the attractor or in the repellor. iv) The gaps in each 
of these families can be ordered with decreasing width in accord to power laws and are seen to 
appear sequentially in the dynamics generated by uniform distributions of initial conditions, v) The 
power law with log-periodic modulation associated to the rate of approach of trajectories towards 
the attractor (and to the repellor) is explained in terms of the progression of gap formation, vi) The 
relationship between the law of rate of convergence to the attractor and the inexhaustible hierarchy 
feature of the preimage structure is elucidated. 
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INTRODUCTION 



In the last few years there has been a reappraisal and 
further exploration of the special dynamical properties 
displayed by critical attractors in simple low-dimensional 
maps This increased attention has been induced in 
part by the perception that the exploration of possible 
limits of validity of the canonical statistical mechanics 
can benefit from the study of much simpler dynamical 
systems that are known to exhibit statistical-mechanical 
analogies 0. For this purpose in mind an ideal model 
system is a one-dimensional map at the transition be- 
tween chaotic and regular behavior, represented by well 
known critical attractors, such as the Feigenbaum attrac- 
tor 3- So far, analytical studies have concentrated 
on the dynamics inside the attractor, characterized by 
fluctuating, memory-preserving, nonmixing, phase-space 
trajectories. These studies have revealed that these tra- 
jectories obey remarkably rich scaling properties. The 
results are exact and clarify [Jl the relationship between 
the original modification 0] [6|_of the thermodynamic 
approach to chaotic attractors [3] [1] ® for this type of 
incipiently chaotic attractor, and some aspects of the q- 
statistical formalism [Tfj| that have been seen to manifest 
for the same critical attractors. Formerly, dynamics in- 
side the attractor was not known at this level of detail 
pXj| . The complementary part of the dynamics, that of 
advance on the way to the attractor, has, to our knowl- 
edge, not been analyzed, nor understood, with similar 
degree of thoroughness. The process of convergence of 
trajectories into the Feigenbaum attractor poses several 
interesting questions that we attempt to answer here and 
elsewhere based on the comprehensive new knowledge 



presented below. Prominent amongst these questions is 
the nature of the connection between the two sets of dy- 
namical properties, within and outside the attractor. 

Trajectories inside the attractor visit positions form- 
ing oscillating deterministic patterns of ever increasing 
amplitude. However, when the trajectories are observed 
only at specified times, positions align according to power 
laws, or g-exponential functions that share the same q- 
index value [4J [lfj. Further, all such sequences of po- 
sitions can be shifted and seen to collapse into a single 
one by a rescaling operation similar to that observed for 
correlations in glassy dynamics, a property known as 'ag- 
ing' [13] (HI- The structure found in the dynamics is 
also seen to consist of a family of Mori's g-phase transi- 
tions via which the connection is made between the 
modified thermodynamic approach and the ^-statistical 
property of the sensitivity to initial conditions Q [13] • 
On the other hand, a foretaste of the nature of the dy- 
namics outside the critical attractor can be appreciated 
by considering the dynamics towards the so-called super- 
cycles, the family of periodic attractors with Lyapunov 
exponents that diverge towards minus infinity. This in- 
finite family of attractors has as accumulation point the 
transition to chaos, which for the period-doubling route 
is the Feigenbaum attractor. As described in Ref. pT| . 
the basins of attraction for the different positions of the 
cycles develop fractal boundaries of increasing complex- 
ity as the period-doubling structure advances towards the 
transition to chaos. The fractal boundaries, formed by 
the preimages of the repellor, display hierarchical struc- 
tures organized according to exponential clusterings that 
manifest in the dynamics as sensitivity to the final state 
and transient chaos. The hierarchical arrangement ex- 
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pands as the period of the supercycle increases [ll| . 

Here we present details on the general procedure fol- 
lowed by trajectories to reach the Feigenbaum attractor, 
and its complementary repellor. We consider an ensem- 
ble of uniformly distributed initial conditions xq spanning 
the entire phase space interval. This is a highly struc- 
tured process encoded in sequences of positions shared by 
as many trajectories with different xq. There is always a 
natural dynamical ordering in the xq as any trajectory of 
length t contains consecutive positions of other trajecto- 
ries of lengths t — 1, t — 2, etc. with initial conditions x' , 
x'q, etc. that are images under repeated map iterations 
of xq. The initial conditions form two sets, dense in each 
other, of preimages of each the attractor and the repellor. 
There is an infinite- level structure within these sets that, 
as we shall see, is reflected by the infinite number of fam- 
ilies of phase-space gaps that complement the multifrac- 
tal layout of both attractor and repellor. These families 
of gaps appear sequentially in the dynamics, beginning 
with the largest and followed by other sets consisting of 
continually increasing elements with decreasing widths. 
The number of gaps in each set of comparable widths 
increases as2 fc , fc = 0, 1, .. . and their widths can be or- 
dered according to power laws of the form a~ k , where 
a is Feigenbaum's universal constant a ~ 2.5091. We 
call k the order of the gap set. Furthermore, by consid- 
ering a fine partition of phase space, we determine the 
overall rate of approach of trajectories towards the at- 
tractor (and to the repellor). This rate is measured by 
the fraction of bins W(t) still occupied by trajectories at 
time t [l2| . The power law with log-periodic modulation 
displayed by W(t) [H[ is explained in terms of the pro- 
gression of gap formation, and its self-similar features are 
seen to originate in the unlimited hierarchy feature of the 
preimage structure. 

Before proceeding to expand our description in the fol- 
lowing sections we recall Q the general definition of the 
interval lengths or diameters d,N, m that measure the bi- 
furcation forks that form the period-doubling cascade se- 
quence in the logistic map f^{x) = 1 — fJ,x 2 , — 1 < x < 1, 
< fi < 2. These quantities are measured when consid- 
ering the superstable periodic orbits of lengths 2 N , i.e. 
the 2 Af -cycles that contain the point x = at JI N < fioo, 
where /ioo — 1.401155189 ... is the value of the control 
arameter /z at the period-doubling accumulation point 
The positions of the limit 2°°-cycle constitute the 
Feigenbaum attractor. The dM.m in these orbits are de- 
fined (here) as the (positive) distances of the elements 
x n , m = 0, 1,2, . . . ,2 N — 1, to their nearest neighbors 



of the Feigenbaum attractor. We present explicit results 
for the logistic map, that has a quadratic maximum, but 
the results are easily extended to unimodal maps with 
general nonlincarity z > 1. 

Central to our discussion is the following broad prop- 
erty: Time evolution at /ioo from t = up to t — > oo 
traces the period-doubling cascade progression from /i = 
up to /ioo. Not only is there a close resemblance be- 
tween the two developments but also quantitative agree- 
ment. For instance, the trajectory inside the Feigenbaum 
attractor with initial condition xq = 0, the 2°°-supercycle 
orbit, takes positions x t such that the distances between 
appropriate pairs of them reproduce the diameters d,N,m 
defined from the supercycle orbits with ~p N < fi^. See 
Fig. 1, where the absolute value of positions and loga- 
rithmic scales arc used to illustrate the equivalence. This 
property has been key to obtain rigorous results for the 
sensitivity to initial conditions for the Feigenbaum at- 
tractor [ij, [3. 
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FIG. 1: Left panel: Absolute value of attractor positions for 
the logistic map / M (x) in logarithmic scale as a function of the 
logarithm of the control parameter difference fioo—f-i- Right 
panel: Absolute value of trajectory positions for the logistic 
map ffj,(x) at /ioo with initial condition xo= in logarithmic 
scale as a function of the logarithm of time t, also show by 
the numbers close to the circles. The arrows indicate the 
equivalence between the diameters djv in the left panel, and 
position differences Dn with respect to xo= in the right 
panel. 
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For large TV, djv,o/<2iV+i,o — a - Notice that infinitely 
many other sequences of superstable attractors appear 
at the period-doubling cascades within the windows of 
periodic attractors for values of /i > /ioo. At each of 
these period-doubling accumulation points occur replicas 



II. PREIMAGE STRUCTURE OF ATTRACTOR 
AND REPELLOR 

A convenient way to visualize how the preimages for 
the Feigenbaum attractor and repellor are distributed 
and organized is to consider the simpler arrangements 



for the preimages of the supercycles' attractors and re- 
pcllors. These arrangements have been characterized in 
detail recently [ll[ for supercycles of small periods 2^, 
N = 1, 2, . . ., and it was observed how they become more 
complicated very rapidly as N grows For N = 1 

the preimages of the attractor belong to only two basins, 
one for each position of the attractor, separated by the 
positions of the fixed-point repellor and its single preim- 
age. For N = 2 the attractor preimages are distributed 
into four basins, one pair of them separated from the 
other pair by a fractal boundary built around the posi- 
tions of the first or 'old' N — 1 repellor and its preimage. 
This fractal boundary consists of the infinite number of 
preimages that the old repellor has acquired and this in 
turn cluster exponentially around it and around its 'old' 
preimage. For N — 3 the attractor preimages are dis- 
tributed into eight basins, four of the basins are sepa- 
rated from the remaining four by a more complex fractal 
boundary sector, built again about the positions of the 
old N — 1 repellor and its old preimage. This sector 
is made of an infinite number of replicas of the fractal 
boundary for the case with N = 2, the increasingly more 
numerous preimages of the old repellor form now an in- 
finite set of clusters, grouping themselves exponentially 
around it and around its old preimage. The basins for 
the positions of the attractor are also separated from each 
other in pairs via other fractal boundary segments with 
the characteristics of the N — 2 case, structures that ap- 
pear centered now around the next generation of repel- 
lors (two) and their first preimages (two). In the panels 
of Fig. [5] we show the relative times t f required to reach 
the attractor as a function of the initial condition xo for 
the supercycles with N = 1, 2, and 3, respectively. The 
logarithmic horizontal scale in these figures relative to 
the position of the old repellor reveals the main charac- 
teristic of the preimage structure for each case. These 
figures illustrate our description above. See Ref. [ll[ for 
the definition of tf and further details. 

As the period 2 N increases the preimage structures 
for the attractor and repellor become more and more in- 
volved, with the appearance of new features made up of 
an infinite repetition of building blocks each of which is 
equivalent to the more dense structures present in the 
previous 2 JV_1 case. In addition all other structures in 
the earlier 2 N ~ 2 , . . . , 2 1 cases are still present. Thus a 
hierarchical organization of preimages is built upon as N 
increases, and the preimage layout for the Feigenbaum 
attractor and repellor is obtained as the limiting form 
of the rank structure of the fractal boundaries between 
attractor basins, that are formed by the preimages of the 
repellor positions for the families of periodic attractors 
symbolized by the supercycles. The fractal boundaries 
consist of sub-basins of preimages for the attractor po- 
sitions separated by preimages of the repellor positions. 
The sizes of these sub-basins decrease while their num- 
bers increase as N — > oo, with the result that in this limit 
their sizes vanish and the preimages of both attractor 
and repellor positions become two sets - with dimension 
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FIG. 2: Time of flight tf(x), the number of iterations 
necessary for a trajectory with initial condition at x to 
reach an attractor position as a function of the distance, 
in logarithmic scale, between x and the repellor position at 
y\. The top panel is for the twice-iterated map fj^(x), 

= 1, yi ~ 0.6180340. The middle panel for f^(x), 
Jl 2 ~ 1.31070264, j/i ~ 0.571663. And the bottom panel 
for f^(x), Tj 3 ~ 1.38154748, yi~ 0.56264475. The values of a; 
near the peaks (open circles) correspond to initial conditions 
very close to the repellor preimages. See text. 

equal to the dimension of phase space - dense in each 
other. (The attractor and repellor are two multifractal 
sets with dimension d f ~ 0.538 ... @|). For the 2 N -th su- 
percycle the positions for its repellor and its first preim- 

(2 N ) 

age appear located at the inflection points of /A (x), 
and it is in the close vicinity of them that the fractal 
boundaries of basins of attraction for the positions of the 
attractor form. To illustrate how the dense sets of preim- 
age structures for the Feigenbaum attractor and repellor 

(2 N ) 

develop we plot in Fig. 3 the absolute value of dfj^ jdx 
for N — 1, 2, . . . , 4 vs. x. The maxima in this curve cor- 

(2 N ) 

respond to the inflection points of (x) at which the 
repellor positions or their first preimages are located. 



III. SEQUENTIAL OPENING OF PHASE 
SPACE GAPS 

One way wherein the preimage structure described 
above manifests in the dynamics is via the successive for- 
mation of phase space gaps that ultimately give rise to 
the attractor and repellor multifractal sets. In order to 
observe explicitly this process we consider an ensamble 
of initial conditions x spread out uniformly across the 
interval — 1 < xo < 1 and keep track of their positions 
at subsequent times. In Figs. 4 to 6 we illustrate the 
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FIG. 3: The absolute value of df£ \x)/dx, for N = 1,2,3 
and 4, in logarithmic scale as a function of x in the interval 
< x < 1. The proliferation of maxima conveys the develop- 
ment of the hierarchical structure of repellor preimages. See 
text. 



outcome for the supercycles of periods 2 2 , 2 3 and 2 4 , re- 
spectively, where we have plotted the time evolution of 
an ensemble composed of 10000 trajectories. In the left 
panel of each figure we show the absolute value of the 
positions \x t \ vs time t, while, for comparison purposes, 
in the right panel we show the absolute value of |a:| both 



4!, '(x) and vs df£ >/dx 



,(2" 



to facilitate identifica- 



tion of the attractor and repellor positions. The labels 
k = 1, 2, 3, . . . indicate the order of the gap set (or equiv- 
alently the order of the repellor generation set [HI). In 
Fig. 4 (with fi = Jl 2 ) one observes a large gap opening 
first that contains the old repellor (k — 0) in its middle 
region and two smaller gaps opening afterward that con- 
tain the two repellors of second generation (k = 1) once 
more around the middle of them. In Fig. 5 (with /i = JI 3 ) 
we initially observe the opening of a primary and the two 
secondary gaps as in the previous ^ = JI 2 case, but subse- 
quently four new smaller gaps open each around the third 
generation of repellor positions (k = 2). In Fig. 6 (with 
/i = /7 4 ) we observe the same development as before, 
however at longer times eight additional and yet smaller 
gaps emerge each around the fourth generation of repel- 
lor positions (k — 3). Naturally, this process continues 
indefinitely as N — > oo and illustrates the property men- 
tioned before for //oo, that time evolution at fixed control 
parameter value resembles progression from ft = up to, 
in this paragraph, Jl N . It is evident in all Figs. 4 to 6 
that the closer the initial conditions Xo are to the repel- 
lor positions the longer times it takes for the resultant 
trajectories to clear the gap regions. This intuitively ev- 
ident feature is essentially linked to the knowledge we 
have gained about the fractal boundaries of the preim- 
age structure, and the observable 'bent over' portions 
of these distinct trajectories in the figures correspond to 



their passage across the boundaries. (Since the ensemble 
used in the numerical experiments is finite there appear 
only a few such trajectories in Figs. 4 to 6). 
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FIG. 4: Phase-space gap formation for /i =~p 2 - Left panel: 
time evolution of a uniform ensemble of 10000 trajectories as 
a function of \x\ (black areas and open circles). The values 
of the index k label the order of the gap set. Right panel: 

Turned around plots of fj^J (a;) (grey) and d/i*^ (x) /da; I (black) 

vs | as | as guides for the identification of attractor and repellor 
positions. 



To facilitate a visual comparison between the process of 
gap formation at /Zoo and the dynamics inside the Feigen- 
baum attractor (as illustrated by the trajectory in Fig. 
lb) we plot in Fig. 7 the time evolution of the same 
ensemble composed of 10000 trajectories with fi = n^. 
We use this time logarithmic scales for both \x t \ and t 
and then superpose on the evolution of the ensemble the 
positions for the trajectory starting at xq = 0. It is clear 
from this figure that the larger gaps that form consec- 
utively have all the same width in the logarithmic scale 
of the plot and therefore their actual widths decrease 
as a power law, the same power law followed, for in- 
stance, by the position subsequence Xt = a~ N , t = 2 , 
N = 0,1,2,... for the trajectory inside the attractor 
starting at xo = 0. This set of gaps develop in time 
beginning with the largest one containing the k = re- 
pellor, then followed by a second gap, one of a set of two 
gaps associated to the k = 1 repellor, next a third gap, 
one gap of a set of four gaps associated to the k = 2 re- 
pellor, and so forth. The locations of this specific family 
of consecutive gaps advance monotonically towards the 
most sparse region of the multifractal attractor located 
at x — 0. The remaining gaps formed at each stage con- 
verge, of course, to locations near other regions of the 
multifractal but are not easily seen in Fig. 7 because of 
the specific way in which this has been plotted (and be- 
cause of the scale used). In Fig. 8 we plot the same data 
differently, with the variable In \x\ replaced by In |1 — x\ 
where now another specific family of gaps, one for each 
value of k = 0, 1, 2, . . ., appear all with the same width in 
the logarithmic scale, their actual widths decrease now as 
a~ 2N , N = 0, 1, 2, . . . The locations of this second family 
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of consecutive gaps advance monotonically towards the 
most crowded region of the multifractal attractor located 
at x = 1 . The time necessary for the formation of succes- 
sive gaps of order k = 0, 1, 2, . . ., increases as 2 fc because 
the duration of equivalent movements of the trajectories 
across the corresponding preimage structures involve the 

2 -th composed function /I ; (x) . 
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FIG. 5: Phase-space gap formation for /i =Jl 3 . Left panel: 
time evolution of a uniform ensemble of 10000 trajectories as 
a function of \x\ (black areas and open circles). The values of 
the index k label the order of the gap set. Right panel: Turned 

around plots of (x) (grey) and \dfj^ (x) /dx I (black) vs \x\ 

as guides for the identification of attractor and repellor posi- 
tions. 
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FIG. 7: Phase-space gap formation for = /loo. The black 
dots correspond to time evolution of a uniform ensemble of 
10000 trajectories as a function of \x\ vs t, both in logarithmic 
scales. The open circles are the positions, labeled by the 
times at which they are reached, for the trajectory inside the 
Feigenbaum attractor with initial condition xo= 0, same as 
right panel in Fig. 1. 
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FIG. 8: Same as Fig. 7 but with replacement of |i| by 
1 — x\. Notice the change in slope with respect to Fig. 7 in 
the opening of gaps and in the layout of the positions for the 
trajectory inside the attractor. 



IV. SCALING FOR THE RATE OF 
CONVERGENCE TO THE ATTRACTOR AND 
REPELLOR 



FIG. 6: Phase-space gap formation for fi = 7J 4 . Left panel: 
time evolution of a uniform ensemble of 10000 trajectories as 
a function of \x\ (black areas and open circles). The values of 
the index k label the order of the gap set. Right panel: Turned 
around plots of /i* 6 ' (a;) (grey) and jd/i.^ 6 - 1 (a;)/da;j (black) vs 

| a; | as guides for the identification of attractor and repellor 
positions. 



There is |12j | an all-inclusive and uncomplicated way to 
measure the rate of convergence of an ensemble of trajec- 
tories to the attractor (and to the repellor) that consists 
of a single time-dependent quantity. A partition of phase 
space is made of Nb equally sized boxes or bins and a uni- 
form distribution, of N c initial conditions placed along 
the interval — 1 < x < 1, is considered again. The num- 
ber r of trajectories per box is r = N c /Nf,. The quantity 
of interest is the number of boxes W(t) that contain tra- 
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jectories at time t. This is shown in Fig. 9 in logarithmic 
scales for the first five supercycles of periods 2 1 to 2 5 
where we can observe the following features: In all cases 
W(t) shows a similar initial nearly constant plateau and 
a final well-defined decay to cero. As it can be observed 
in the left panel of Fig. [9] the duration of the final de- 
cay grows (approximately) proportional to the period 2 N 
of the supercycle. There is an intermediate slow decay 
of W{t) that develops as N increases with duration also 
(just about) proportional to 2 N . For the shortest period 
2 1 there is no intermediate feature in W(t), this appears 
first for period 2 2 as a single dip and expands with one 
undulation every time N increases by one unit. The ex- 
panding intermediate regime exhibits the development 
of a power-law decay with the logarithmic oscillations 
characteristic of discrete scale invariance [lij. Clearly, 
the manifestation of discrete invariance is expected to be 
associated to the period-doubling cascade. In the right 
panel of Fig. [5]we show a superposition of the five curves 
in Fig. [5] (left panel) obtained via rescaling of both W(t) 
and t for each curve according to repeated scale factors. 

The limiting form W(t) for N — > oo is shown in the 
left panel of Fig. [TU] for various values of r while in its 
right panel we show, for r = 100, a scale amplification of 
W{t) with the same factors employed in Fig. [9] for the 
supercycles with small periods. The behavior of Wit) at 
fioc was originally presented in Ref. [12j where the power 
law exponent s and the logarithmic oscillation parameter 
A in 

Wm.H^y- ,2) 

were obtained numerically with a precision that corre- 
sponds to r — 10. In Eq. (2) h(x) is a periodic function 
and A is the scaling factor between the periods of two con- 
secutive oscillations. More recently, in Ref. [l5| it was 
pointed out that numerical estimates of W{t) are sub- 
ject to large finite-size corrections, and, also, that W(t) 
should scale with the intervals in the triadic cantor set 
construction of the Feigenbaum attractor [l5[ [l6| , from 
which the value for s = 0.800138194 is reported. The 
values for the rescaling factors in our Figs. M and [TUJ suf- 
fer from these large finite size effects due to the relatively 
small values of r used in the calculations. This is evident 
since the time scaling factor obtained from these data 
differs 10% from the exact value of A = 2 implied by the 
discrete scale invariance property created by the period- 
doubling cascade. In Fig. [TJJwe show the rate W(i) and 
the superposition of repeated amplifications of itself (as 
in the right panel of Fig. [TO)) for increasing values of N c . 
We find that the scaling factor A converges to its limit 
A = 2. 

We are now in a position to appreciate the dynamical 
mechanism at work behind the features of the decay rate 
W(t). From our previous discussion we know that every 
time the period of a supercycle increases from 2 JV_1 to 2^ 
by shifting the control parameter value from /Zjv-i to ~p N , 
the preimage structure advances one stage of complica- 




FIG. 9: Left panel: The rate W(t), divided by the number 
of boxes Nb employed, of approach to the attractor for the 
supercycles of periods 2 , N = 1, 2, 3,4 and 5 in logarithmic 
scales. The expression shown corresponds to the power-law 
decay of the developing logarithmic oscillations. Right panel: 
Superposition of the five curves for Wit) in the left panel via 
n-times repeated rescaling factors shown for the horizontal x 
and vertical y axis. 
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FIG. 10: Left panel: The rate W(i) of approach to the Feigen- 
baum attractor at fi = /i^ in logarithmic scales. The curves 
correspond to the values given for the number r of trayec- 
tories per box, and the expression shown corresponds to the 
power-law decay with logarithmic oscillations. Right panel: 
Superposition of Wit), for r = 100, with itself via the rescal- 
ing shown (the same as in Fig. 8) for the horizontal x and 
vertical y axis. 



tion in their hierarchy. Along with this, and in relation to 
the time evolution of the ensemble of trajectories, an ad- 
ditional set of 2 N smaller phase-space gaps develops and 
also a further oscillation takes place in the corresponding 
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rate Wit) for finite-period attractors. At fJ, = /loo time 
evolution tracks the period-doubling cascade progression 
and the flow of trajectories undergo every time t increases 
from 2 N ~ 1 to 2 N equivalent passages across stages in 
their itinerary through the preimage ladder structure, in 
the development of phase-space gaps, and in logarithmic 
oscillations in W(t). In Fig. 12 we show the correspon- 
dence between the latter features quantitatively. Also, in 
this figure we have added the results of a calculation of 
W(t) at times t = 2 N , N = 0, 1, 2, . . ., according to the 
expression 

N 

W(2 N ) = ]T d N , m . (3) 

m— 1 

This confirms the value s = 0.8001. 
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FIG. 11: The same as in the right panel of Fig. 10 but 
obtained with increased precision. In the left panel the num- 
ber of initial conditions is N c = 10 6 , while in the right panel 
7V C = 10 7 . The distance A is A = (1 + | - l/a|)/(l + | - 1|), 
where a is Feigenbaum's constant. This stems from the fact 
that all initial conditions out of the interval ( — 1/a, 1) take a 
value inside this interval in the first iteration. As it can be 
observed the scaling factor for the horizontal axis converges 
to the exact value x — 2. 



V. SUMMARY 

We have examined the process followed by an ensem- 
ble of uniformly distributed initial conditions xq across 
the phase space to arrive at the Feigenbaum attractor, or 
get captured by its corresponding repellor. Significantly, 
we have gained understanding concerning the dynamical 
ordering in the xq , in relation to the construction of the 
families of phase-space gaps that support the attractor 
and repellor, and about the rate of approach of trajec- 
tories towards these multifractal sets, as measured by 




FIG. 12: Correspondence between the power-law decay with 
log-periodic oscillation features of the rate W(t) with the se- 
quential opening of phase space gaps. Top panel: The solid 
line is Wit) from Fig. 9 and the open circles values are ob- 
tained for Wit) from Eq. (3) at times t = 2 N , N = 1,2,... 
See text. 



the fraction of bins W(t) still occupied by trajectories at 
time t. An important factor in obtaining this knowledge 
has been the consideration of the equivalent dynamical 
properties for the supercycles of small periods in the bi- 
furcation cascade [ll| . 

As we have seen, a doubling of the period introduces 
well-defined additional elements in the hierarchy of the 
preimage structure, in the family of phase space gaps, 
and in the log-periodic power law decay of the rate W(t). 
We have then corroborated the wide-ranging correla- 
tion between time evolution at /ioo from t = up to 
t —>■ oo with the 'static' period-doubling cascade pro- 
gression from \i = up to \ioo. As a result of this we 
have acquired an objective insight into the complex dy- 
namical phenomena that fix the decay rate Wit). Wc 
have clarified the genuine mechanism by means of which 
the discrete scale invariance implied by the log-periodic 
property in Wit) arises, that is, we have seen how its self- 
similarity originates in the infinite hierarchy formed by 
the preimage structure of the attractor and repellor. The 
rate W(t) can be obtained quantitatively (see Eq. (3)) 
from the supercycle diameters d^.m- This basic data de- 
scriptive of the period-doubling route to chaos is also a 
sufficient ingredient in the determination of the anoma- 
lous sensitivity to initial conditions for the dynamics in- 
side the Feigenbaum attractor [Ioj |. 
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